Plasmons and near-field amplification in double-layer graphene 
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We study the optical properties of double-layer graphene for linearly polarized evanescent modes 
and discuss the in-phase and out-of-phase plasmon modes for both, longitudinal and transverse 
polarization. We find a energy for which reflection is zero, leading to exponentially amplified trans- 
mitted modes similar to what happens in left-handed materials. For layers with equal densities 
n = 10 12 cm -2 , we find a typical layer separation of d ~ 500/im to detect this amplification for 
transverse polarization which may serve as an indirect observation of transverse plasmons. When 
the two graphene layers lie on different chemical potentials, the exponential amplification either fol- 
lows the in-phase or out-of-phase plasmon mode depending on the order of the low- and high-density 
layer. This opens up the possibility of a tunable near-field amplifier or switch. 

PACS numbers: 73.21. Ac, 42.25.Bs, 78.67.Wj, 73.20.Mf 
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I. INTRODUCTION 

The optical properties of graphene have attracted 
much attention recently. 1 A hallmark is given by the uni- 
versal absorption of 2.3% in the visible regime indepen- 
dent of the incident wavelength. 2,3 For infrared frequen- 
cies, the absorption can be altered by varying the Fermi 
energy via an applied gate voltage 4 . For fixed chemical 
potential, the optical conductivity is then characterized 
by three regimes: i) the Drude peak related to intraband 
transitions, ii) a region with suppressed absorption due to 
Pauli-blocking and iii) the nearly constant plateau with 
oo = § tt due to interband transitions. 5,6 

There is also high potential to use graphene in 
information and communication technologies. Apart 
from applications as conductive and optically trans- 
parent electrode in resistive touchscreens of transpar- 
ent and flexible displays, 7,8 there are various exam- 
ples to use graphene for opto-electronic devices such as 
photodetector 10 , broadband absorber 9 , or mode- locked 
laser due to current saturation 11 . There are also inter- 
esting theoretical proposals not yet observed like the dy- 
namical generation of a gap 12 or carrier multiplication 
due to strong electron-electron interaction 13 . 

Here, we analyze the optical properties of double- 
layer graphene which is of potential interest due to the 
recent advances in fabricating graphene-boron nitride 
structures. 14,15 Our main result is that there is exponen- 
tial amplification of evanescent modes in gated double- 
layer systems of graphene, first discussed by Pendry in 
the context of so-called left-handed materials with neg- 
ative refraction index. 16 The appearance of exponential 
amplification in between two layers is closely related to 
the existence of surface plasmon polaritons as noted by 
Pendry and later discussed by Haldane 17 . This effect can 
be considered the evanescent equivalent of the extraor- 
dinary transmission observed in perforated metal sheets, 
where plasmons are believed to be involved. 18-20 

In graphene, there are two types of (surface) plasmons, 
longitudinal (or transverse magnetic) 21,22 and transverse 



(or transverse electric) modes 23 which give rise to dif- 
ferent behavior. Theoretically, these low-energy collec- 
tive excitations have attracted much interest 24-31 and 
recently, it was shown that they can be excited via a 
attenuated total reflection structure, 32 paving the way 
to graphene plasmonics. 33-35 Experimentally, they have 
been observed via electron energy loss spectroscopy 36,37 
and near-field nanoscopy. 38 

Here, we will discuss how the plasmon modes are modi- 
fied in a double-layer structure for both polarizations. In 
the case of longitudinal polarization, the two plasmon 
modes will split, if the layer separation is small enough 
such that the plasmon modes can still considerably in- 
teract electrostatically. This leads to an ordinary two- 
dimensional plasmon with Jq dispersion but larger en- 
ergy where the charges of the two layers oscillate in phase, 
and a linear (acoustic) plasmon mode where the charges 
oscillate out of phase. Their dependence on the Fermi 
energy of the two layers has first been discussed in Ref. 39 
for T = 0, for finite temperature in Ref. 40 , and the plas- 
mon interaction and hybridization in pairs of neighboring 
nanoribbons in Ref. 41 . Here, we also discuss the scat- 
tering problem. In the case of transverse polarization, 
plasmons are not formed via charge accumulation and 
are only weakly bound to the graphene layer. These col- 
lective excitations are linked to the two-band structure 
of the electronic carriers and thus also exist in bilayer 
graphene where they become more localized. 42 In a dou- 
ble layer structure, we find either one mode (for small 
layer separation) or two plasmon modes (for large layer 
separation) . 

As noted before, in the case of two plasmon modes, 
there exists a frequency oj exp where reflection is zero and 
the transmission is exponentially amplified. For equal 
electronic densities in the two layers, this frequency lies 
between the two plasmon frequencies and is pinned to the 
out-of-phase mode for small layer separation. For large 
layer separation, the two plasmon modes merge and sym- 
metrically sandwich the frequency of exponential amplifi- 
cation just as in the original model proposed by Pendry. 16 
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For different electronic densities, Lo exp can follow either 
the low- or high-density plasmon mode depending on the 
order of the layers. 

We also discuss the amplification A proportional to 
the ratio of the intensity at the second interface to one 
at the first interface. The maximum of A lies between the 
two plasmon frequencies but only coincides with uj exp in 
the case of large layer separation. This tunable plasmonic 
evanescent amplification could be useful for near-field mi- 
croscopies. 

The paper is organized as follows. In section II, we 
present the scattering problem and discuss the general 
behavior. In section III, we discuss the case of equal 
electronic density in the two layers and in section IV, we 
discuss the density-imbalanced system. We close with a 
summary and conclusions. In an appendix, the explicit 
formulas for the general setup with different dielectric 
media are given. To model graphene, we will only de- 
pend on the current-current correlation function at zero 
temperature which was evaluated in Ref. 44 within the 
Dirac cone approximation and in Ref. 45 for the hexago- 
nal lattice. 



II. SCATTERING PROBLEM 

To discuss the scattering behavior of linearly polarized 
light in a layered geometry, a convenient gauge condition 
is to choose the scalar potential to be zero. We thus only 
have to consider the vector potential. 

In the following, we will assume the graphene layers to 
be parallel to the zy-plane. For transverse polarized light 
with say q = qe x , only the y-component A y is non-zero. 
Since here we are interested in the scattering behavior 
of evanescent modes, we defin e the real transverse wave 
number q[ = \/ q 2 — tiiuj j c) 2 , where q denotes the con- 
served wavenumber parallel to the plane and the di- 
electric constant of the corresponding medium. For two 
dielectric media separated by graphene located at z = 0, 
the vector potential of an incoming field reads 
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A similar expression holds for longitudinal polarization, 43 
see also the appendix. 

For transverse polarization, the vector potential is con- 
tinuous at the interface and the first derivative makes a 
jump due to the current generated by the vector field in- 
side the graphene plane. A similar situation occurs for 
longitudinal polarization (see appendix). The current is 
related to the corresponding current-current susceptibil- 
ity ^H'- 1 -, 44 ' 45 where || stands for longitudinal or p and 
_L for transverse or s polarization. This results in the 



following expressions for the transmission amplitude: 

+ ^ + 

T \\ = 2ei^ 

ei<7 2 + £ 29i 7-p 



(2) 
(3) 



where £o is the electric permittivity of free space. We 
further assumed the same magnetic permittivity for the 
two regions. The reflected mode is given by R = T — 1. 

The general case of several interfaces can be obtained 
from the above result by summing up the (infinite) tra- 
jectories leading to the total transmission and reflection, 
respectively. The explicit expressions for two graphene 
layers sandwiched by three dielectric media are presented 
in the appendix, solving the matching conditions directly. 



A. Homogeneous medium 

In the following, we will discuss the special case of a 
homogeneous medium where the three dielectrica are the 
same in all three regions e = ej. This allows us to treat 
both polarizations on the same footing. The general case 
is discussed in the appendix. 

For two graphene layers separated by the distance d 
and using the notation of Eq. (1), we obtain 

1 



R 



1 + ai + a 2 + aia 2 (l - e- 2q ' d ) 1 
ai + a 2 e- 2( i' d + a x a 2 (\ - e - 2q ' d ) 



(4) 
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l + a 1 +a 2 +a 1 a 2 {l-e- 2c i' d ) ' 

where at — xi / {^ e oc 2 q') for transverse and at = 
~x!V/(2e£o^ 2 ) f or longitudinal polarization. 



At the energy fuu exp with 

ot\ + a 2 e~ 2q d + aia 2 (l 



e- 2q ' d ) = , 
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there is no reflection, R — 0, the evanescent equivalent to 
perfect transmission. At this energy, we have exponential 
amplification of the transmission amplitude 



T - - — 

± exp 



(7) 

a 2 a 2 

where we introduced the (scalar) conductivities Uj of the 
two layers well defined for Fermi energies less than leV. 45 
We will also discuss the amplification related to the ra- 
tio of the two intensities at the two interfaces, see Eq. 
(14). The maximum is then centered around R = — 1 
where the expression divergences. This corresponds to 
the condition 

l + a 2 (l-e- 2 i' d ) 



. 



(8) 

Let us now introduce the dielectric matrix e(q, u) 
which for a homogeneous medium is given by 39 



1 + ai ct\e 



-q d 



a 2 e- q ' d l + a 2 
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Genuine plasmons correspond to the zeros of the deter- 
minant dete, 

e(g,w) = (1 + a x )(l + a 2 ) - a 1 a 2 e~ 2q ' d . (10) 

This quantity was discussed in Ref. 39. For two un- 
damped plasmon modes (lmai=0), the slopes around the 
two zeros are opposite and the corresponding delta func- 
tions obtained by hne~ 1 (q, to + iO) thus have weights of 
different signs. In fact, also in the dissipative region the 
two modes can be related to opposite signs of the so- 
called "loss" function lme~ 1 (q,u> + iO), contrary to the 
results presented in Ref. 39. 

The sign change of e(q,ui) docs not violate causality 
since it is only the determinant of a 2 x 2-matrix response 
function. Apart from its zeros, it does not have any fur- 
ther physical interpretation. In this work, we thus prefer 
to discuss the eigenvalues of the 2 x 2-matrix current- 
current response function defined as j = -e 2 xA el ' where 
j = (ji , j 2 ) T denotes the current in layer one and two. We 
obtain the following expression: 

— / \ -1 ( ai(l + a 2 ) aia 2 e~ q ' d \ n , 

X(9,w =— ,J , (11 

de{q,Lo) \ a x a 2 e qa a 2 {l + ot 1 )J 

where d = q' e 2 / (2c£quj 2 ) for longitudinal polarization 
and d — —e 2 /{2e a c 2 q') for transverse polarization. In 
the following, we will discuss the sum of the two eigen- 
values or, equivalently, the trace of x(g,w), 

X(q,u) = , i \ («i + ot-2 + 2a x a 2 ) . (12) 
ae{q,u)) 

We will call — Imx(q, w + iO) the energy loss function as 
it provides the total plasmonic spectral weight. 

III. LAYERS WITH EQUAL DENSITIES 

Let us discuss some basic properties of the double-layer 
system in a homogeneous medium for the case of equal 
electronic densities, i.e., a = a\ = a 2 . This situation is 
easily realized by applying a potential difference between 
the two layers which results in an equal density of elec- 
trons and holes in the two layers, respectively. Since the 
carrier type is irrelevant for the optical properties as long 
as particle-hole symmetry is retained, this experimental 
setup is equivalent to the one with equal polar doping. 

The plasmonic excitations are given by the zeros of 
e(q, ui), see Eq. (10). For lma=0, this yields the following 
conditions satisfied by the plasmon frequency ui^(q), 

l + a(l±e-«' d ) = 0, (13) 

where the plus (minus) sign refers to the in-phase (out- 
of-phase) plasmon mode. We will use these conditions 
to find resonances also in the case of so-called Landau 
damping (Ima 7^ due to particle- hole excitations) , thus 
formally setting Ima = also in the region of interband 
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FIG. 1: (color online): The in-phase (upper blue curve) and 
out-of-phase (lower blue curve) longitudinal plasmon mode 
for two different layer separations d — 2, lOnm are shown 
together with the single layer plasmon mode obtained for 
d — > 00 (dashed-dotted curve). Also show are the energies 
corresponding to R = (dashed magenta curve) and R = — 1 
(dashed red curve). Density of the two graphene layers is 

ml2 —2 j.x-_ J: r„._ „ 
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FIG. 2: (color online): The in-phase (upper blue curve) and 
out-of-phase (lower blue curve) transverse plasmon mode for 
two different layer separations d — 50, 100/im are shown 
together with the single layer plasmon mode obtained for 
d — y 00 (dashed-dotted curve). Also show are the energies 
corresponding to R = (dashed magenta curve) and R = — 1 
(dashed red curve). Density of the two graphene layers is 
ni — n-z = 10 12 cm -2 with the same dielectric medium for the 
three regimes e = 1 (vacuum) and we set Ima = 0. 



transitions (lo/vf > 2kp — q). For large layer separation 
d and/or large wave number q' ', the two plasmon modes 
thus converge to the plasmon condition for a single layer. 
This is expected from physical grounds since for q'd 3> 1, 
the two layers are decoupled and each layer has the same 
plasmonic properties. In Refs. 39,40 the imaginary part 
of a is considered throughout and the two branches dis- 
perse when they cross the region of interband transitions. 
Both methods (setting lma=0 or not) produce similar re- 
sults for weak damping, the only relevant limit to still be 
able to speak of a collective excitation. Still, we have 
found from comparison with — Imx(g,w + iO) that set- 
ting Ima = also in the region of interband transitions 
provides a better procedure to obtain the maximum of 
the plasmonic spectral function. Also the energies cor- 
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responding to R — and R = — 1 will be obtained by 
setting Inw = 0. 

For q'd 3> 1, Eq. (6) coincides with the plasmon con- 
dition for uncoupled layers (d — ¥ oo), i.e., 1 + a = 0. For 
q'd <C 1 , we have for the out-of-phase mode the simplified 
condition 1 + aq'd = 0. In this limit, this coincides with 
the condition for zero reflectivity given in Eq. (6). The 
branch of uj exp thus shifts from to w+ as the layer 
separation d is increased. For d — > oo, uj exv is located in 
between the two plasmonic resonances ojp just as in the 
original model of Pendry with a left-handed lens. 46 The 
reason why this is not the case for small layer separation 
is just because the two oscillators have different oscillator 
strengths which shifts the resonance condition. 

In Fig. I (longitudinal polarization) and Fig. 2 (trans- 
verse polarization), the in-phase (upper blue curve) and 
out-of-phase (lower blue curve) plasmon mode for two 
different layer separations d are shown together with the 
single layer plasmon mode obtained for d — > oo (dashed- 
dotted curve). We also show the energies correspond- 
ing to R = and R = —I. For both polarizations, 
we choose n\ = = 10 12 cm -2 which corresponds to 
kp O.I77nm -1 and the same dielectric medium for the 
three regimes with e = I (vacuum). 

As already discussed, in the above figures we have only 
considered the real part of the susceptibility and thus ne- 
glected the possibility of damping. At zero temperature, 
damping will occur at wave numbers q > 2kp — u)/vp 
due to interband transitions. This is indicated by one 
of the dotted lines in the two figures. In Fig. 1, also 
the limiting dispersion for the acoustic mode u> = VFq is 
shown. Even though the "approximate" analytical for- 
mula of Ref. 39 for longitudinal polarization states that 
the slope is proportional to the Fermi wave vector which 
could be chosen arbitrarily small, the square-root singu- 
larity of the susceptibility at u> — vf q forces the acoustic 
plasmon velocity to remain greater than vp. 

For transverse polarization, usually the long wave- 
length approximation of the susceptibility is used as- 
suming an infinite speed of light. This yields a loga- 
rithmic divergence for uj/vpkp = 2. 23 Including the full 
g-dependence, we see two finite maxima at w/uf&f = 
2/(1 ± vf/c) when the light cone crosses the regions de- 
fined by the straight lines lo/vf = 2k f ± q- 

A. Damping and amplification of evanescent modes 

To discuss the damping of the plasmonic mode, we plot 
the energy loss function, — Imx(g r , a; + iO). The results 
are shown on the left hand side of Fig. 3 (longitudinal 
polarization) and Fig. 4 (transverse polarization) . Intro- 
ducing a infinitesimal imaginary part to the frequency, 
uj — s- w + iO, gives rise to delta functions defining the plas- 
monic frequencies in the region where there is no damp- 
ing. These are shown as solid blue lines, which broaden 
for q > 2kp — us/vp. 

Let us now discuss the amplification of evanescent 




FIG. 3: (color online): Left hand side: The energy loss func- 
tion — Imx{q,(^ + iO) in units of h = Vp = 1 for longitudinal 
polarization. Right hand side: The amplification A of lon- 
gitudinal polarization together with the plasmonic dispersion 
of the in-phase and out-of-phase plasmon mode (solid blue 
curves). Also the energies corresponding to R — (dashed 
magenta curve) and R = — 1 (dashed red curve) are shown. 
The density of the two graphene layers is n\ = ri2 = 10 12 cm -2 
with the same dielectric medium for the three regimes e — 1 
(vacuum). 



modes. For this, we introduce A proportional to the 
ratio of the field intensity at the second interface with 
respect to the one at the first interface: 



The amplification diverges for R = — 1 and the corre- 
spding energy lies in between the two plasmon modes 
as seen from Eq. (8), but does not coincide with the 
exponential amplification given in Eq. (6) (only in the 
limit d — ¥ oo, the two expressions converge). This can 
be seen on the right hand side of Figs. 3 (longitudinal 
polarization) and 4 (transverse polarization), where the 
amplification A is shown together with the energies corre- 
sponding to R — (dashed magenta curve) and R = — 1 
(dashed red curve). We note that even after including a 
small finite damping term, the amplification is strongly 
peaked around the red dashed curve (R = —1) and only 
the grafical cutoff makes A appear to be smeared out. 

For large layer separation, though, the two plasmon 
modes coincide and also the maximum of A is centered 
around this curve as well as uj exp . This is not shown here 
in the case of equal densities, but the situation is similar 
to the one shown on the right hand side of Fig. 6 and 
Fig. 7. 

B. Detecting transverse plasmons 

For transverse polarization, the plasmon modes are 
close to the light cone. 23 We, therefore, introduced the 
adimensional variable Q = cq/(vpkp) and showed the 
difference Q — iojvpkp in Fig. 2 and 4. The fact that 
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FIG. 4: (color online): Left hand side: The energy loss func- 
tion —lm.x(q,oj + iO) in units of % = Dp = 1 for transverse 
polarization. Right hand side: The amplification A for trans- 
verse polarization together with the plasmonic dispersion 
of the in-phase and out-of-phase plasmon mode (solid blue 
curves). Also the energies corresponding to R = (dashed 
magenta curve) and R = — 1 (dashed red curve) are shown. 
The density of the two graphene layers is n\ — ri2 = 10 12 cm -2 
with the same dielectric medium for the three regimes e = 1 
(vacuum) . 



they are almost pinned to the light cone makes them dim- 
cult to observe. Recently, it was proposed to use strained 
graphene to detect them, since strain slightly changes the 
dispersion relation. 48 An alternative detection method is 
based on fluorescence quenching of graphene which is en- 
tirely determined by transverse plasmons in an appropri- 
ate energy window. 49 Here, we want to stress that also 
the double-layer geometry could be used to indirectly de- 
tect these collective excitations which are linked to the 
two-band structure of the carriers of graphene. 

The main idea is that exponential amplification is 
linked to the two plasmon modes. By increasing the 
interlayer distance d, the amplification can become ar- 
bitrarily large, see Eq. (7). But for large separation, the 
interaction between the two layers decreases and dissipa- 
tive effects will eventually render the effect. 

Let us conclude with some numbers. For n = 
10 12 cm~ 2 , the layer separation should be larger than 
d>50^m, the distance where the second (out-of-phase) 
mode appears. For larger separation d>lmm, the two 
plasmon modes almost coincide, and the exponential am- 
plification w exp lies almost on top of these modes. We 
would thus expect an optimal separation to lie around 
the two limiting separations, i.e., d ~ 100 ~ 500/im, 
suitable to observe strong amplification. 



IV. LAYERS WITH DIFFERENT DENSITIES 

Let us now assume a density difference between the 
two layers. For longitudinal plasmons, we will investi- 
gate the case where one layer has a typical density of 
n\ = 10 12 cm -2 and the other layer a (slightly) lower 
doping level with ri2 — 10 11 cm -2 . We note that even 



at the neutrality point there are (damped) plasmons at 
finite temperature which can be treated with the ana- 
lytical formulas given above for finite kp — > 7^2 In 2. 47 
At room temperature, this corresponds approximately to 
n = 10 cm -2 , but in this case also strong damping terms 
are present. 

Due to the lesser lateral confinement of the plasmon 
on the layer with lower density, larger layer separations 
are now possible until convergence of the two plasmon 
modes sets in. By convergence we now mean that the 
two plasmon modes are defined by the free conditions 
1 + oti =0 and 1 + 012 — in their respective energy 
window without intcrband transitions. In Fig. 5, we 
show the energy loss function for two different layer sep- 
arations d — 2nm and d — 50nm. There are thus two 
regimes: in the first regime with d<20nm, there are an 
acoustic and a plasmonic-like mode. Both modes are 
only undamped in the region where no interband transi- 
tions are present. When interband transitions due to the 
lower density layer set in, both modes are damped but 
the plasmonic one reaches up to the region where inter- 
band transitions due to the higher density set in. For the 
second regime with <i>20nm, there is a almost free plas- 
mon with 1 + ct\ rs for q > 2k p — uj/vf limited only by 
interband transitions related to the higher density. Also 
the free plasmon 1 + a>i — is shown as dashed black line 
as guide for the eye. One the other hand, the acoustic 
mode has become plasmonic with 1 + a.2 ~ 0, limited by 
interband transitions related to the lower density. 

The condition for exponential amplification depends 
on whether the electromagnetic wave first passes the high 
density layer or the low density layer, see Eq. 6. In Fig. 
6, the amplification A is shown for the first case, choosing 
the layer densities as ri\ — 10 12 cm~ 2 and n% — 10 cm -2 . 
The amplification A follows the out-of-phase mode and 
terminates at the region where interband transitions re- 
lated with the lower density set in. 

Also shown is the condition for exponential amplifica- 
tion (magenta curve for Imcti = 0) which yields a low 
and a high-energy solution. The low energy solution 
follows the low energy plasmon mode for both regimes 
(left and right panel). The high energy solution starts 
at uj/{vpkp) = 1.667 and then disperses towards higher 
energy roughly following the plasmonic in-phase mode. 

In Fig. 7, the the amplification A together with the 
plasmon modes (setting Ima^ = 0) is shown for the 
opposite layer densities, thus choosing ni = 10 11 cm -2 
and ri2 = 10 12 cm -2 . For (very) small layer separation 
d <C 20nm, the maximum of A is centralized in the vicin- 
ity of the out-of-phase mode (not shown), then crosses 
the region between the two plasmonic modes (left hand 
side), until it reaches the in-phase mode for large layer 
separation d>20nm (right hand side). 

The amplification now only terminates at the region 
defined by interband transitions related with the higher 
density. This can be understood from Eq. (8), which only 
contains the physical properties of the second (higher- 
energy) layer. This condition is shown for Imcv2 = as 
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q/k F ' ' q/k F 

FIG. 5: (color online): Left hand side: The energy loss func- 
tion — Imx(g,^ + iO) in units of h = Vp = 1 for longitudinal 
polarization and two different layer separations d = 2nm (left) 
and d = 50nm (right). Right hand side: Also shown the plas- 
monic dispersion of the in-phase plasmon mode as obtained 
without the imaginary part of Qi (black). The density of the 
two graphene layers is n\ — 10 12 cm -2 and ri2 = 10 11 cm -2 
with the same dielectric medium for the three regimes e = 1 
(vacuum) . 




FIG. 7: (color online): The amplification A for longitudi- 
nal polarization and two different distances d = 2nm (left) 
and d = 200nm (right) choosing n\ = 10 11 cm -2 and ri2 = 
10 12 cm -2 . The plasmonic modes (black) conditions for ex- 
ponential amplification (magenta) with Ima?i = are also 
shown. Additionally, on the left hand side, we present the 
energy corresponding to R — — I (dashed red curve). 




FIG. 6: (color online): The amplification A for longitudi- 
nal polarization and two different distances d — 2nm (left) 
and d — 200/^m (right) choosing ri2 = 10 11 cm -2 and m = 
10 12 cm -2 . The plasmonic modes (black) and the conditions 
for exponential amplification (magenta) with Ima^ = are 
also shown. 



dashed red curve on the left hand side of Fig. 7 and is 
slightly shifted towards higher energies than the maxi- 
mum of A due to the neglect of damping terms. Again 
we stress that A has a pronounced resonance close to 
the energy corresponding to R = —1 and only appears 
smeard out due to the chosen grafical cutoff. 

The condition for exponential amplification (magenta 
curve for Imaj = 0) now yields different solutions de- 
pending on the layer separation. For d<20nm, we again 
find a low and a high-energy solution. For the second 
regime d>20nm, the two solutions for exponential am- 
plification split. There is a branch which only exists at 
small wavevectors and within an energy window reaching 
from uj/(vpkp) ~ 1.667 to ui/(vpkp) 2. The second 
branch starts at w/(vpkp) « 2 for q = and then fol- 



lows the in-phase mode when it crosses the region of in- 
terband transitions coming from the lower density layer. 
The third branch first follows the in-phase mode for small 
q wave vectors and then disperses towards smaller en- 
ergies at the (avoided) crossing point with the second 
branch. We note that this crossing happens at the en- 
ergy ui/(vFkp) « 1.667. The repeated occurrence of the 
value u/(vFk F ) « 1.667 is due to the fact, that at this 
energy the response of the ith graphene sheet becomes 
zero, which obviously affects the response of the other 
graphene sheet. 

Let us close with a comment on transverse plasmons. 
For this polarization, there is generally only a small 
window in which undamped plasmons can exist, i.e., 
1.667 < O < 2 with £1 = w/(w F fcF)- 23 Since for expo- 
nential amplification, the two plasmons of the two layers 
have to exist at the same frequency, this limits the density 
difference to kp/kp<1.2. But from Fig. 2, we see that 
the plasmon with lower energy exists in an even smaller 
energy window and we obtain solutions for typical densi- 
ties n = 10 12 cm -2 only for kp/kp < 1.02. Nevertheless, 
there is a rich phenomenology also in the case of a larger 
density difference. We leave this analysis when detailed 
experimental setups are available. 



V. SUMMARY AND CONCLUSIONS 

We have discussed the optical properties of double 
layer graphene for linearly polarized evanescent modes 
and focused on the plasmon dispersion and (exponential) 
amplification of evanescent modes. We considered longi- 
tudinal as well as transverse polarization and especially 
the detection of exponential amplification for transverse 
polarized modes might at last yield evidence of trans- 
verse plasmons, uniquely related to the chirality of the 
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FIG. 8: (color online): Schematic setup of the double-layer 
graphene structure. The two graphene layers, characterized 
by ai and a-z and separated by a distance d, are sandwiched 
by the three dielectric media characterized by ei, £2, and £3. 
The income wave (normalized to one) is reflected (R) and 
transmitted (T). 



electronic carriers of graphene and not found in other 
materials. 

For systems with different densities for the two layers, 
we observe that the transmission crucially depends on the 
the order of the two layers. If the first layer is the high 
density layer, the amplification of the evanescent modes 
follows the dispersion of the low energy solution, whereas 
for the other case, the high energy solution is relevant. 
By varying the relative density of the two layers, one 
can switch from amplifications based on the out-of-phase 
mode to ones based on the in-phase mode. This might 
give rise to interesting applications where the modulation 
of an external gate can be used for an optical switcher 
or, more generally, for a tunable near field amplifier. 
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Appendix A: General formulas for three different 
dielectrica 



Here, we present the general formulas of the reflectivity 
and transmissivity for a structure with three different 
dielectrica and graphene in there interfaces, i.e., there 
are three dielectric media separated by graphene located 
at z — and z — d. The graphene layers shall be parallel 
to the xy-plane. A schematic representation of the setup 
is shown in Fig. 8. We will discuss the scattering problem 
only for the vector potential. We thus work with a gauge 
in which the scalar potential is set to zero. 



For transverse polarized light with say q = qe x , only 
the y-component A y is non-zero. Since we are interested 
in the time-evolution of evanescent modes, we define the 
real transverse wave number q[ = y 'q 2 — ei(uj/c) 2 , where 
q denotes the conserved wavenumber parallel to the plane 
and ei dielectric constant of the corresponding medium. 



1. Transverse polarization 

The vector potential for a transversed polarized filed 
with q = qe x is given by 



A y (r,z)=J2e iqxX 



e -9i« 4 

ae~^ z 
Te- q 's z 



Re q i z 
f be q 'i z 



,z < 

, < z < d 

,z> d 



(Al) 

The vector field is continuous at the interface. From 
V x B = /ioj, it follows that the first derivative of A y 
is discontinuous at the interface due to the current of the 
form j = jie y S(z — Zj) with z^ = 0,d. Since the current 
is generated by the vector field, i.e., fa — —\^A y {zi) 
where xt is the current-current response function of the 
graphene layer at z = Zi, the set of equations closes and 
we obtain the following expressions: 



2g[(q' 2 + q' 3 + 2a 2 ) 
N 

2e-^ d q' 1 (q' 2 -q' 3 ~2a 2 ) 
N 



(A2) 

(A3) 
(A4) 



N = (q[ +q 2 + 2 ai ){q' 3 + q' 2 + 2a 2 ) 

-( q ' 1 -q 2 +2a 1 )(q 3 -q 2 + 2a 2 )e-^ d 

where a; = x^/(2soc 2 ). This leads to 

T 1 - = Aq' iq ' 2 e {q '^ )d /N , R ± = a + b-l . (A5) 

At the energy hoj exp with 

(q , 2 -q[ + 2a 1 ){q' 3 + q , 2 + 2a 2 ) (A6) 

=(q' 1 +q' 2 -2a 1 )(q 2 -q 3 -2a 2 )e- 2q * d , 

there is no reflection, B = 0. At this energy, we have 
exponential amplification of the transmission amplitude 



, q 2 -q' 1 + 2a 1 (q > +q > )d 
"** q' 2 -q' 3 ~2a 2 e 
The amplification A diverges for 

(q' 2 +q' 3 + 2a 2 ) = -(q' 2 -q' 3 - 2a 2 )e~ 2q '- d 
2. Longitudinal polarization 



(A7) 



(A8) 



For longitudinal polarization, the general vector field 
also has a component in normal direction to the inter- 
faces: 

A(r, z) = J2 eiq ' r (^11 (q> ^) e q + ^±(q. 2 K) (A9) 
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where with q[ = yj q 2 — €i(uj/c) 2 we make the ansatz 

(*=IU) 



e- q i z + R l e q ^ z 
a^r^ + bie q '^ 



,z<0 

, < z < d 

,z>d 



(A10) 



The components A± are obtained from A\\ via the con- 
dition for a transverse field V • A = 0. In the following 
we will drop the index i and evaluate the coefficients for 

The vector field is continuous at the interface, but 
the normal component of the displacement field makes 
a jump due to the presence of graphene and the corre- 
sponding charge density in the layer. This component is 
related to the vector field via the relation D z = e eiuiA z . 
With the continuity equation up — q • j = and the 
linear response j = — xf-Aq, the set of equations closes. 
Assuming the same notation as in the case of transverse 
radiation, we obtain the following expressions: 

_ 25261(9362 + q' 2 e 3 + 2a 2 ) 



N 



b = 



2e-^ d q 2 e 1 (q' 3 e 2 - q ' 2 e 3 -2a 2 ) 



(All) 
(A12) 

N = {qi.d + q[e 2 + 2a 1 ){q' 2 e 3 + q' 3 e 2 + 2a 2 ) (A13) 
- (<7 2 ei - + 2 ai ){q' 2 e 3 - q' 3 e 2 + 2a 2 )e- 2q 'i d 



N 



where ce\ = — x{ q[q 2 / (2eou 2 ) and a 2 = 
-X [ lq 2q l 3 /{2e uj 2 ). This leads to 

Til = iq'^e^e^'^/N , R.W = a + b - 1 . (A14) 

At the energy haj exp with 

(q[e 2 - q 2 t x + 2a x ){c( 2 e z + q' 3 e 2 + 2a 2 ) (A15) 
= (<Z 2 ei + q [e 2 - 2 ai )(q' 3 e 2 - q' 2 e 3 - 2a 2 ) e - 2 ^ d , 



there is no reflection, R = 0. At this energy, we have 
exponential amplification of the transmission amplitude 



T ll _ q[e 2 - g 2 ei + goi q' 3 {q > 2+q > 3)d 
£XP &2 - q' 2 e 3 - 2a 2 q[ 



(A16) 



The amplification A diverges for 

(93 e 2 + q' 2 €3 + 2a 2 ) = -(q' 3 e 2 - q 2 e 3 - 2a 2 )e~ 2q2<i . 

(A17) 
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